Statistical Mechanics of segregation in binary granular mixtures 
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We discuss mixing/segregation phenomena in a schematic hard spheres lattice model for binary 
mixtures of granular media, by analytical evaluation, within Bethe-Peierls approximation, of Ed- 
wards' partition function. The presence of fluid-crystal phase transitions in the system drives segre- 
gation as a form of phase separation. Within a pure phase, gravity can also induce a kind of vertical 
segregation not associated to phase transitions. 
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Granular mixtures subject to shaking can mix or, un- 
der different conditions, spontaneously segregate their 
components according to criteria which, although of deep 
practical and conceptual importance, are still largely un- 
clear 0,0. In this respect, microscopical grains proper- 
ties are known to play a role since, for example, smaller 
grains appear to filter beneath larger ones (the so called 
"percolation" phenomena 0, 0). "Inertia" 0, "con- 
vection" and "dynamical" effects in general are also 
known to be relevant [jj . Interestingly, recent simulations 
and experiments have outlined that segregation processes 
can also involve "global" mechanisms: phenomena such 
as "condensation" [8j or, more generally, "phase separa- 
tion" s| have been found, recalling well known proper- 
ties in the Statistical Mechanics of thermal fluids, such 
as colloidal mixtures (see Ref.s in 0,0]). Actually, even 
though granular media are not characterized by Boltz- 
mann distributions fluidized mixtures have been 

schematically treated as standard fluids under gravity in 
presence of a finite bath temperature 0] . 

A different Statistical Mechanics approach was intro- 
duced by Edwards 0,0,0 (see for a review). This 
was explicitly designed to deal with non-thermal systems, 
such as granular media in their "jammed states" (i.e., not 
in their fluidized regime): time averages of a granular 
system (subject to some drive, e.g., tapping) observed 
at rest, are supposed to coincide with suitable ensem- 
ble averages over its "mechanically stable" states, i.e., 
those where the system is found still. Even though the 
limit of validity of such an approach must be still as- 
sessed, it appears to be well grounded in several cases 
0,0,0,0]. In particular, it was shown to hold in a 
schematic model of granular mixtures: a lattice binary 
hard spheres system under gravity subject to sequences 
of taps ^ij. In this paper we analytically solve, at the 
level of Bethe approximation, the partition function a la 
Edwards of such a model and derive its phase diagram 
and mixing/segregation properties as a function of grains 
masses, sizes, number and, see below, "configurational 
temperatures" . 

Our model exhibits phase transitions from fluid to crys- 
tal phases. As much as in thermal media, this induces 
segregation effects associated to the presence of phase 
separation phenomena, with the formation of coexisting 



phases rich in small or large grains. Gravity drives ver- 
tical segregation, where large grains can be found on av- 
erage above (the well known "Brazil nut" effect, BNE 
ll|) or below small grains (reverse BNE, RBNE 0,0). 
When associated to phase separation, we have a form 
of strong segregation, as opposed to mixing, or similarly 
weak segregation, observed within a given phase. Coars- 
ening phenomena, as those experimentally observed 

0, 

are predicted to occur for suitable values of the system 
parameters. 

The model under investigation |l4| is a binary hard 
sphere mixture made of two species, 1 (small) and 2 
(large) with grain diameters qq = 1 and V2ao, under 
gravity on a cubic lattice, of spacing <zo, confined in a 
rigid box. On each lattice site we define an occupancy 
variable: n\ — 0, 1, 2 if site i at height z is empty, filled by 
a small or by a large grain. The system Hamiltonian is: 
H = Hhc + migHx + m 2 gH 2 , where Hx = J2i, z zS nfi, 
H 2 = J2i z z &n*2 are the heights of the two species and 
TLhc is the hard core potential, preventing two nearest 
neighbor sites to be both occupied if at least one contains 
a large grain. 

As stated, Edwards' approach appears to hold with 
good approximation in this model j 1 11] under a tap dy- 
namics (analogous to that used in the compaction of 
real granular media), which turns out to be character- 
ized by two suitable "thermodynamic parameters" , con- 
jugated to the heights of the two species and called 
configurational temperatures, T^ on ^ and T£° . Their 
values are increasing functions of the tap amplitude, 



tap 

r, and tend to coincide for r — > |l5j. In brief, 
the weight of a given microstate r turns out to be 
0: exp{-H H c(r) ~ PwngH^r) - fom 2 gH 2 {r)} • IF, 



where T^ onf = /?f 1 and T. 
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(3 2 ■ The operator IF r 



selects mechanically stable states: II r = 1 if r is "stable" , 
else IF = 0. The system partition function a la Edwards 
is thus the following |l4|: 

2 = ^""^ e -[HHc(r)+l3im 1 gH l +li 2 m-2gH-2\ . pj^ 

M 

where the sum is over all microstates r. Since the exact 
calculation of Z is hardly feasible, we want to evaluate 
Z here at a mean field level. To this aim we consider 
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a generalization of Bethe-Peierls method for anisotropic 
systems (due to gravity), i.e., we solve by recurrence re- 
lations the partition function of the system on the Bethe 
lattice shown in the inset of Fig^ (see also 0, El ) . The 
reader not interested in technicalities about calculations 
can jump to the paragraph after Eq.Q where the results 
are discussed. In particular, we consider a 3D lattice 
box with H horizontal layers (i.e., z <E {1, ...,H}) occu- 
pied by hard spheres. Each layer is a random graph of 
given connectivity, k — 1 (we take k = 5). Each site in 
layer z is also connected to its homologous site in z — 1 
and z + 1 (the total connectivity is thus k + 1). The 
Hamiltonian is the one above plus two chemical potential 
terms to control the two species densities (//i,/^)- Hard 
Core repulsion prevents two particles on connected sites 
to overlap. In the present lattice model we adopt a simple 
definition of "mechanical stability": a grain is "stable" 
if it has a grain underneath. For a given grains config- 
uration r = {n,}, the operator II r has a tractable ex- 
pression: IL, = liniR-^oo exp{-KHEdw} where Tisdw = 



0n !: 9U Z — l n O z-2 

n i z n. On. 



The local tree-like properties of our lattice allows to 
write down recursive equations a la Bethe found by iter- 
ation of the lattice structure. Details on calculations and 
notations can be found in 0,0,0], here we just recall 
the main lines. A "branch" of our lattice can iteratively 
grow in the "up", "down" or "side" directions. The par- 
tition functions of new branches starting from a site i at 
height z can be recursively written in terms of their old 

(i z) (i z) 

ones. To this aim define Zq ' s and Zn's the partition 
functions of the "side" branch restricted respectively to 
configurations in which the site i is empty or filled by a 
particle of specie n (with n 



7 {i ' z) and 7 (M 
^0,« ana A n,u 



1,2); analogously, Z%£', 



(resp. Z ( ^\ Z { Q l f and Z^l) are the 
partition functions of the "up" (resp. "down") branch 
restricted to configurations in which the site i is filled by 
a grain of specie n, empty with the upper (resp. lower) 
site empty and empty with the upper (resp. lower) site 
filled by a grain of specie n. In order to write Bethe re- 
cursive equations, it is convenient to introduce the local 
"cavity fields" defined by: e s ™ 



,(«,*) 



(»,*) /y (i,z) „(v) 



In these new vari- 
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e^' 2) =^Sl4f (with n= 1,2) 
ables the recursion relations are more easily written 
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Here n = 1,2, the products run on the neighbors, j, of 
= l + On,ie 1 , T± J ' = 1 + e"i + e ^ , 
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FIG. 1: Main frame Phase diagram, from Bethe approxi- 
mation, of the hard sphere binary mixture model confined to 
a single 2D horizontal layer. Here the control parameters are 
the densities per unit surface of small and large grains Ni, 
N2. A Fluid and a Crystal phase are found, divided by a 
discontinuous melting transition. When the mixture is pre- 
pared with values of iVi and N2 out of the Fluid or Crystal 
phases (i.e., in the white area), it exhibits phase separation, 
associated to coarsening, in a fluid phase rich in small grains 
and a crystalline phase rich in large ones. This causes hori- 
zontal segregation of the two grains species. We also plot the 
metastable "supercooled" fluid region, where a coexistence 
line between two supercooled fluids, with different relative 
species concentrations, is found (dashed line) having a critical 
point (open dot). For high N2, a 1-step Replica Symmetry 
Breaking glassy phase (described in details elsewhere [l9|p. 
metastable with respect to the crystal, is present. Inset A 
schematic picture of the 3D lattice model: we consider a bi- 
nary mixture of hard spheres located on a Bethe lattice where 
each horizontal layer is a random graph of given connectivity. 
Homologous sites on neighboring layers are also linked and 
the overall vertex connectivity isc = fc + l = 6. 

From local fields the free energy, F, can be derived 
H> El along w hh the quantities of interest, such 
as the density profile of small and large grains, p\(z) 
and P2(z), their density per unit surface Ni = p\{z) 
and N2 = X] z P2(z) and average heights h n = (z n ) = 
J2z z Pn( z )/J2 z Pn( z ) (with n = 1,2). The system pa- 
rameters (for a given grains sizes ratio) are four: the two 
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number densities per unit surface, iVi and N2 (conju- 
gated to the chemical potentials), and the two configu- 
rational temperatures, or more precisely m\(5\ and to 2 /32 
(conjugated to gravitational energies). In the space of 
these parameters, the fluid phase corresponds to a solu- 
tion of Bethe-Peierls equations where local fields in each 
layer are site independent. Such a solution, character- 
ized by horizontal translational invariance, is given by the 
fixed points of Eqs.J2J. With the standard Bethe-Peierls 
methods, a crystalline phase, characterized by the break- 
down of the translational invariance (local fields are now 
different on neighboring sites), is also found. The typical 
resulting phase diagrams are shown in Fig.sn an d|3 
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FIG. 2: Phase diagram of a 3D system under gravity, treated 
a la Edwards, in the plane (Ni,N2) for mi/3i = 1 and 
mifii = 1/2 (main frame). A pure Fluid and Crystal phase 
are present along with a region where a form of phase separa- 
tion occurs and the two phases coexist (the white area marked 
Crystal&Fluid): since gravity breaks the vertical symmetry, 
here a crystal rich in large grains is resting on a fluid bed 
rich in small grains. This is a strong BNE situation, visu- 
alized in the right inset showing the density profiles of the 
two species, pi(z) and £2(2) (resp. filled and empty circles), 
in a typical point of the Crystal&Fluid region with Ni — 3, 
Na = 4. The reverse, i.e., strong RBNE with the fluid float- 
ing above the crystal, can be found when mi/3i < ma^i- 
When the configurational temperatures and species densities 
get close coarsening phenomena may occur. For comparison 
in the left inset we show p\(z) and pi(z) as a function of 
the vertical coordinate 2, in a point of the Crystal phase with 
Ni = 0.3, N2 = 4: small grains are here interspersed with 
large ones even though, on average, slightly below. This illus- 
trates that within a pure phase, gravity also drives a weaker 
form of segregation, not associated to phase separation, as 
discussed in the text. For clarity, metastable phases are not 
shown in this phase diagram. 

For sake of clarity, we first describe the case where the 
system is confined to a 2D horizontal layer (see Fig. QJ. 
In this case gravity plays no role (i.e., we can set g = 0) 
and the only control parameters are Ni and -/V2. Our 



findings are consistent with known results on binary hard 
spheres with no gravity [2(J. A Fluid phase is found for 
small N2 and, by increasing N 2 , there is a discontinuous 
(except for iVi = 0) melting transition to a Crystal phase. 
Whenever the system is prepared with values of N\ and 
N2 out of the regions where these phases are stable (i.e., 
in the white area of Fig|TJ, phase separation occurs be- 
tween a fluid rich in small grains and a crystal rich in 
large grains |2l|. This might correspond to a horizontal 
segregation in the grains mixture, as recently found in 
the form of coarsening phenomena in experiments on a 
monolayer by Reis and Mullin 0. Fig. ^ a lso shows the 
"supercooled fluid" and "glassy" metastable phases 
(plotted with smaller fonts). Since nucleation times can 
be in practice very long (and enhanced by a degree of 
polidispersity), crystallization can be avoided and the 
metastable fluid observed. In such a region we find a 
coexistence line between a fluid rich in small grains and 
another one rich in big grains (dashed line in Fig. ^| with 
a second order critical point (big empty circle). 

Fig. |21 shows the typical phase diagram of a 3D sys- 
tem in a case where mifti > m 2 P2- The pure Fluid and 
Crystal phases found in 2D are still present, even though 
their extension depends now on /3± and (3% (they shrink as 
/3i and P2 increase) . There is also a region where the two 
phases coexist (marked Crystal&Fluid). Gravity breaks 
the system "up-down" symmetry, so the coexisting sta- 
ble or metastable phases of the 2D case are vertically 
segregated. This separation mechanism underlies a form 
of strong vertical segregation (see Fig|2J): here the crystal 
phase, rich in large grains, moves to the top and a clear 
cut BNE is found (as RBNE is observed in the oppo- 
site case, when m\{i\ < m^fa)- This kind of segregation 
is visualized by a plot of the density profiles which, in 
this region, show a clear separation of the two coexisting 
phases (see right inset of Fig 0). 

Opposed to such phase separation driven segregation, 
within the pure Fluid and Crystal phases one observes 
mixing or a form of weak vertical segregation since more 
mechanically stable states can be found with small grains 
at the bottom (as described by "percolation" mecha- 
nisms). This is shown, in a typical point of the Crys- 
tal phase, by the species density profiles plotted in the 
left inset of FigEl small grains are essentially mixed 
with large ones and, on average, slightly below. In the 
phase diagram of Fig |3 different forms of BNE can thus 
be found in different regions (similarly for RBNE in the 
case mi Pi < miP-i). In general, for a given grains sizes 
ratio, an interplay of mass densities and configurational 
temperatures difference drives the phases vertical posi- 
tioning. 

Due to the symmetry breaking field, in 3D no coars- 
ening phenomena are usually associated to segregation, 
in contrast with the 2D case. Coarsening is expected to 
appear in the segregation process when both the phases 
densities and the configurational temperatures get close, 
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a phenomenon which could be confirmed by experiments 
or simulations. 

These are the basic mechanisms which underly seg- 
regation/mixing phenomena and which give rise to a 
variety of behaviors as grains and system parameters 
are changed. This is illustrated by a few more exam- 
ples in Fig|21 selected to describe cases where the sys- 
tem moves from weak BNE to weak RBNE regions. The 
picture shows the usual vertical segregation parameter 
Ah/h = 2{h\ — h 2 )/(hi + h 2 ). In the main panel of Fig0 
Ah/h is plotted as a function of N 2 in a region where the 
system crosses from the fluid to the crystal phase (see 
Fig0), the full dot showing the phase transition point: 
well within the crystal, Ah/h changes its sign passing 
from BNE, Ah/h < 0, to RBNE, Ah/h > 0. The insets 
of FigOHshow Ah/h in the fluid phase as a function of 
Ni and T™ nf = T 2 zonf = T con f respectively. 
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FIG. 3: In a 3D system a variety of mixing/segregation be- 
haviors can be observed by changing grains and system pa- 
rameters. Main frame Given N\ = 0.4, mi/3i = 0.8 and 
wi2,02 = 1.25, the species relative height difference Ah/h is 
plotted as a function of N2 as the system crosses from the 
Fluid to the Crystal phase (see Fig|5J, the full dot and the 
dashed line showing the phase transition point. Upper inset 
Given N2 = 0.4, TOi/?i = 1 and m2/32 = 2, Ah/h is plotted 
as a function of Ni in the Fluid phase. Lower inset Given 
Ni = 0.5, N 2 = Ni/2, mi = 1 and m 2 = 2, Ah/h is plotted 
as a function of T 1 con/ = T 2 con/ = T conf in the Fluid phase. 
In all these cases Ah/h changes sign, passing from a BNE to 
a RBNE region, through a mixing region where Ah/h ~ 0. 

Summarizing, the present mean field Statistical Me- 
chanics model of granular mixture, here analytically 
treated a la Edwards, allows to explain the basic mech- 
anisms underlying a variety of mixing and segregation 
phenomena experimentally observed, ranging from BNE 
Hi and RBNE HG3, to 

coarsening effects Si. Interest- 
ingly for non-thermal media, they turn out to be related 
to thermodynamic-like mechanisms taking place in the 
system phase diagram. 



Work supported by MURST-PRIN 2002, MIUR-FIRB 
2002, CRdC-AMRA, INFM-PCI. 



[1] J.M. Ottino and D.V. Khakhar, Ann. Rev. Fluid Mech. 
32, 55 (2000). T. Shinbrot and F.J. Muzzio, Physics To- 
day March issue, 25 (2000). J. Bridgewater, Chem. Eng. 
Sci. 50, 4081 (1995). 

[2] "Unifying concepts in granular media and glasses", (Else- 
vier Amsterdam, in press), Edt.s A. Coniglio, A. Fierro, 
H.J. Herrmann, M. Nicodemi. 

[3] A. Rosato, K.J. Strandburg, F. Prinz, R.H. Swendsen, 
Phys. Rev. Lett. 58, 1038 (1987). 

[4] J. Bridgewater, Powder Technol. 15, 215 (1976). J.C. 
Williams, Powder Technol. 15, 245 (1976). 

[5] T. Shinbrot and F.J. Muzzio, Phys. Rev. Lett. 81, 4365 
(1998). 

[6] J.B. Knight, H.M. Jaeger, SR. Nagel, Phys. Rev. Lett. 

70, 3728 (1993). 
[7] J.T. Jenkins and D.K. Yoon, Phys. Rev. Lett. 88, 194301 

(2002) . 

[8] D.C. Hong, PV. Quinn, S. Luding, Phys. Rev. Lett. 86, 
3423 (2001). 

[9] K. M. Hill and J. Kakalios, Phys. Rev. E 49, R3610 
(1994). P.M. Reis and T. Mullin, Phys. Rev. Lett. 89, 
244301 (2002). P.M. Reis, G. Ehrhardt, A. Stephenson 
and T. Mullin, Europhys. Lett. 66, 357 (2004). S. Au- 
maitre, T. Schnautz, C.A. Kruelle, and I. Rehberg, Phys. 
Rev. Lett. 90, 114302 (2003). 
[10] M. Fasolo and P. Sollich, Phys. Rev. Lett. 91, 068301 

(2003) . 

[11] A. Kudrolli, M. Wolpert, J. P. Gollub, Phys. Rev. Lett. 
78, 1383 (1997). J.S. Olafsen and J.S. Urbach, Phys. Rev. 
Lett. 81, 4369 (1998). F. Rouyer and N. Menon, Phys. 
Rev. Lett. 85, 3676 (2000). G.W. Baxter and J.S. Olaf- 
sen, Nature 425, 680 (2003). 

[12] J. A. Both and D.C. Hong, Phys. Rev. Lett. 88, 124301 
(2002). 

[13] S.F. Edwards and R.B.S. Oakeshott, Physica A 157, 
1080 (1989). A. Mehta and S.F. Edwards, Physica A 157, 
1091 (1989). 

[14] M. Nicodemi, A. Fierro, A. Coniglio, Europhys. Lett. 60, 
684 (2002). 

[15] A. Fierro, M. Nicodemi, A. Coniglio, J. Phys.: Condens. 

Matter 15, S1095 (2003). A. Fierro, M. Nicodemi, A. 

Coniglio, to be published. 
[16] M. Nicodemi, Phys. Rev. Lett. 82, 3734 (1999). A. Barrat 

et al, Phys. Rev. Lett. 85, 5034 (2000). J.J. Brey, A. 

Prados, B. Sanchez-Rey, Physica A 275, 310 (2000). A. 

Coniglio and M. Nicodemi, Physica A 296, 451 (2001). D. 

S. Dean and A. Lefevre, Phys. Rev. Lett. 86, 5639 (2001). 

H. A. Makse and J. Kurchan, Nature 415, 614 (2002). 

A. Fierro, M. Nicodemi and A. Coniglio, Europhys. Lett. 

59, 642 (2002). G. De Smedt, C. Godreche, J.M. Luck, 

Eur. Phys. J. B 32, 215-225 (2003). 
[17] A. P. J. Breu, H.-M. Ensner, C. A. Kruelle, and I. Re- 
hberg, Phys. Rev. Lett. 90, 014302 (2003) 
[18] M. Mezard and G. Parisi, Eur. Phys. J. B 20, 217 (2001). 

G. Biroli and M. Mezard, Phys. Rev. Lett. 88, 025501 

(2002). 

[19] A. Coniglio, A. de Candia, A. Fierro, M. Nicodemi, M. 



5 



Tarzia, Europhys. Lett. 66 in press; cond-mat/0309612 
[20] see for instance A. Imhof and J.K.G. Dhont Phys. Rev. 
Lett. 75, 1662 (1995). M. Dijkstra, R. van Roij, R. Evans, 
Phys. Rev. E 59, 5744 (1999), and ref.s therein. 



[21] For sake of clarity here we considered the case where 
small grains cannot form a crystalline phase. The latter 
can be obtained with a slightly more complex TLhc- 



